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We present a detailed analytical and numerical study of nonequilibrium 
dynamics for the complex Ginzburg-Landau (CGL) equation. In particular, 
we characterize evolution morphologies using spiral defects. This paper (re- 
ferred to as I) is the first in a two-stage exposition. Here, we present analyti- 
cal results for the correlation function arising from a single-spiral morphology. 
We also critically examine the utility of the Gaussian auxiliary field (GAF) 
ansatz in characterizing a multi-spiral morphology. In the next paper of this 
exposition (referred to as II), we will present detailed numerical results. 
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I. INTRODUCTION 



There have been many studies of pattern formation in the complex Ginzburg-Landau 
(CGL) equation, which has the general form: 

= mt) + (1 +^«)vv(r,t) - {i+zpmf,t)\'mt). a) 

In Eq. (H), ip{r,t) is a complex order-parameter field which depends on space (r*) and 
time (t); and a, (3 are real parameters. The CGL equation arises in diverse contexts, e.g., 
chemical oscillations thermal convection in binary fluids @; multi-mode lasers f^; etc. 
An overview of applications of the CGL equation is provided in the review article by Cross 
and Hohenberg . The importance of the CGL equation stems from the fact that it provides 
a generic description of the slow modulation of oscillations in a spatially-extended system 
near a Hopf bifurcation . 

The CGL equation exhibits a rich range of dynamical behavior with variation of the 
parameters a and (3, and the "phase diagram" has been investigated (mostly numerically) 
by various authors [^]. In a large range of parameter space, the emergence and interaction 
of spiral (and antispiral) defects plays an important role in determining the morphology. 
Our present work focuses on characterizing pattern formation in the CGL equation using 
spiral-defect structures. 

We have analytically and numerically studied nonequilibrium dynamics in the CGL 
equation with a = 0. Typically, we consider the evolution morphology resulting from a 
small-amplitude random initial condition. There has been intense research interest in such 
problems in the context of far-from-equilibrium statistical physics - for reviews, see 
The simplest problem in this class considers a homogeneous two-phase mixture, which has 
been rendered thermodynamically unstable by a rapid quench below the critical co-existence 
temperature. An example of such a system is a 2-state ferromagnet (in zero magnetic 
field) at high temperatures, which consists of a homogeneous mixture of "up" and "down" 
spins. However, below the critical temperature, the system prefers to be in a spontaneously- 
magnetized state. The evolution of the system from the unstable initial state is a complex 
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nonlinear process. In appropriate dimensionless units, this evolution is described by the 
time-dependent Ginzburg-Landau (TDGL) equation, i.e., Eq. (|T]) with ip{r,t) real, and 
a = (3 = 0. The system evolves by the formation and growth of domains which are enriched 
in either up or down spins, and are characterized by a time- dependent length scale, L{t). In 
the case of a pure and isotropic ferromagnet, the domain growth law is L(t) ~ t^/^, which is 
referred to as the Lifshitz-Cahn-Allen (LCA) law The primary mechanism for domain 
coarsening (or "phase ordering dynamics") is the curvature-driven motion and annihilation 
of interfaces (or defects). Ohta et al. [|l^ have formulated an interface-dynamics approach 
to obtain an analytic form for the equal-time correlation function of a phase-ordering ferro- 
magnet. 

Next, let us consider the dynamical XY model, which is Eq. (|l|) with ip{f,t) complex, 
but a = (3 = 0. In this case, the relevant defect structures (for dimensionality d > 2) are 
vortices (or vortex lines, etc.), and domain growth is driven by the motion and annihilation 
of vortices and anti- vortices. Puri has obtained the time- dependent correlation function 
for the XY model, using singular-perturbation methods due to Suzuki [|T2|, Kawasaki et al. 
TB| , and Puri and Roland Furthermore, Bray and Puri |T5| and (independently) Toyoki 



have obtained the time-dependent correlation function for the vector TDGL equation 
with 0{n) symmetry in d dimensions when n < d, i.e., when topological defects are present. 
(The dynamical XY model corresponds to the case with 0(2) symmetry.) The corresponding 
domain growth law is again the LCA law, L{t) ~ t^/^, with logarithmic corrections when 
n = d PJI?!). To the best of our knowledge, there are no general results available for the 
case with n > d, where the absence of topological defects makes it difficult to characterize 
the dynamical evolution. 

The present two-stage exposition focuses on phase ordering dynamics in the CGL equa- 
tion with a = 0. Furthermore, the analytical and numerical results presented here are 
for the 2-dimensional case, where spirals are point defects. However, the analytical results 
obtained by us can be easily extended to the case with a 7^ and c? > 2, as the underly- 
ing paradigm remains the same, i.e., spiral defects still determine the morphology in large 



regions of parameter space and for higher dimensionahty. 

Following the work of Hagan |T^, Aranson et al. |T^, and Chate and Manneville 0, 
we briefly discuss the phase diagram of the d = 2 CGL equation with a = 0. The limit 
(3 = corresponds to the dynamical XY model, which is well understood [p!T],p|. Without 
loss of generality, we consider the case with /? > 0. For < /5 < /?i ~ 1.397 fT^)? spirals 
(which are asymptotically plane-waves) are linearly stable to fluctuations. For I3i < (3 < (32 
(/?2 — 1.82 spirals are linearly unstable to fluctuations, but the growing fluctuations 

are advected away, i.e., the spiral structure is globally stable. Finally, for (32 < (3, the 
spirals are globally unstable structures and cannot exist for extended times . Our results 
correspond to the parameter regime with (3 < (32- 

In this paper (referred to as I), we present analytical studies of the correlation function 
resulting from single-spiral and multi-spiral morphologies. The next paper (referred to as 
II) presents detailed numerical results and compares them with the analytical results of 
I. This paper is organized as follows. In Section II, we obtain analytical results for the 
correlation function of a single-spiral morphology. In Section III, we critically examine the 
utility of the Gaussian auxiliary field (GAF) ansatz for the characterization of a multi- 
spiral morphology. Section IV concludes this paper with a brief summary and discussion of 
our analytical results. 



II. CORRELATION FUNCTION FOR A SINGLE-SPIRAL MORPHOLOGY 

Figure 1 shows a typical evolution from a small-amplitude random initial condition for 
the d = 2 CGL equation with a = and (3 = 1. We have plotted constant-phase regions in 
this figure, and it is clear that the evolving morphology is characterized by spirals and their 
interactions. (We use the term "spiral" for both spirals and antispirals, unless specifically 
stated otherwise.) There is a characteristic length scale, e.g., inter-spiral spacing or square 
root of inverse defect density, which we denote as L. Details of our simulation techniques 
and comprehensive numerical results will be provided in paper II. Figure 1 is shown here 
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only to motivate our subsequent discussion. 

We would like to quantitatively characterize the evolution morphology shown in Figure 1. 
The standard tool for this is the correlation function of the order-parameter field [j^fl, which 
we will define shortly. (The momentum-space structure factor is obtained as the Fourier 
transform of the real-space correlation function.) At the simplest level of approximation, 
the morphology in the frames of Figure 1 can be interpreted as consisting of disjoint spirals, 
each of size L. (Of course, this overlooks modulations of the order-parameter field at spiral- 
spiral boundaries, but we will discuss that later.) Therefore, it is obviously of relevance to 
compute the correlation function for a single-spiral solution. 

The CGL equation with a = has been studied by Hagan [0, who found that there is 
a family of spiral solutions with the following functional form (in d = 2): 

tjj{f,t)=p{r)exp —il3{l — (f)t + im6 — i(j){r) , (2) 

where r = (r, 0); g > is a constant which is determined uniquely as a function of /?; and 
m is the number of arms in the spiral. The cases with m > and m < correspond to a 
spiral and antispiral, respectively. The limiting forms of the functions p(r) and 0(r) are 

p(r) (1 — g^)^^^, (p'ir) — >■ g, as r — >• oo, 

p(r) ar'^, 4>'{r) ^ r, as r — > 0, (3) 

where the constant a is determined by finiteness conditions. Hagan has presented explicit 
solutions for q{(3) in the cases with m = 1, 2. We will focus on the case with m = ±1, as only 
the 1-armed spirals are expected to be stable in the evolution |jl8|. Figure 2 plots Hagan's 



solution for q{(3) (with f3 < 1.5) in the case with m = ±1. In the simple limit /? = 0, we 
have q = 0, and the spiral solution simplifies to the vortex solution - for the m = ±1 vortex, 
the lines of constant phase correspond to constant 6. Spiral solutions for the general case 
with q;,/5 7^ have been discussed by Aranson et al. [ p^ , pO| ] . 

We are interested in the correlation function for a 1-armed spiral at large length scales, 
so we simplify Eq. (||) as 



'0(r, t) ~ yjl — q'^ exp —if3{l — q'^)t + — qr) 



(4) 



where we have speciahzed to the m = 1 case. The correlation between points n and r2 is 
determined as 

C{n,r2,t) = Re{ij{n,t)r{r2,t)} 

~ (1 - g^)Re {exp[i{6i - qri) - i{62 - qr2)]} 

^C{fuf2). (5) 

The average correlation function is obtained by integrating over the point fi, setting r2 = 
ri + fi2, i.e., 



C{ri2) = J driC{fi,fi + ri2)h{L - |ri + fi2|) 



V 



Re / driexp[i{9i - 02 - qri + q\ri + ri2\)]h{L - \ri + rul), 



(6) 



where V is the spiral volume. In Eq. (j^), we use the step function, h{x) = 1 (0) if x > 
(x < 0), which ensures that we do not include points which lie outside the defect of size L. 

For d = 2, the vector notation ^2 = fi + fi2 is equivalent to r2e*^^ = rie*^^ + ri2e*^^^. 
Thus, we have 

rie'^i + ri2e'^i2 



je2 



[rf + ri2^ + 2riri2 cos(6'i - 6*12)] 



1/2 ■ 



(7) 



and 



C(ri2) 



exp ( — < ri 



/•Z/ /'27r 

Re / driri / (i6' 
JO Jo 



ri + rue 



i{ei-ei2) 



[rf + ri2^ + 2riri2 cos(^i - 612)] 



1/2 



1/2 



rf + + 2riri2 cos(^i - ^12) U(L - |fi + fi2|). (8) 



We introduce the variables 9i — 612 = 0; x = ri/L; r = ri2/ L, to obtain 

X + re 



C(ri2) 



^ ^Re / dxx I dd- 



+ + 2xr cos ^)^/2 



exp 



-iqL |x - (x^ + + 2xr cos 6')^/^| h[l - (x^ + + 2xr cos6')^/^] 



(9) 



where we have used V = vrL^ in = 2. Thus, the scaling form of the single-spiral correlation 
function is C(ri2)/C(0) = g{ri2/L,q'^L'^). In general, there is no scaling with the spiral size 
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because of the additional factor qL. We recover scaling only in the limit q = {P = 0), 
which corresponds to the case of a vortex. Essentially, spirals of different sizes are not 
morphologically equivalent because there is more rotation in the phase as one goes out 
further from the core. 

Figure 3 plots C(ri2)/C(0) vs. ri2/L for the case with /3 = 1 (g ^ 0.306). These results 
are obtained by a direct numerical integration of Eq. (P). We consider 4 different values 
of L. The functional form in Figure 3 exhibits near-monotonic behavior for small values 
of L (i.e., in the vortex limit); and pronounced oscillatory behavior for larger value of L, 
as is expected from the integral expression. Notice that ri2/L < 2 - larger values of ri2 
correspond to the point r2 lying outside the defect. 

Before we proceed, we should point out that the imaginary part of the integral in Eq. 
is non-zero, in general - corresponding to a weak correlation between the real and imaginary 
parts of the order-parameter field. The imaginary part can also be obtained with relative 
ease. However, we will confine our discussion to the conventional definition of the correlation 
function in Eq. (^. Let us next consider the asymptotic behavior of the correlation function 
in the limit ri2/L 0, though ri2 is still much larger than the size of the defect core ^. 

A. Case with /? = 

In the case with /3 = 0, we have g = and the integral expression in Eq. @ simplifies 

as 

(~'[ri2) = —Re / dxx / dU-r-^ ttttt 

^ ' vr io Jo (x2 + r2 + 2xrcos^)V2 

The behavior in the r — > limit is of considerable interest as it determines the large- 
wavevector {k — * oo) behavior of the structure factor In that case, we can neglect the 
step function on the right-hand-side (RHS) of Eq. ( p!OD as it only provides corrections at the 
edge of the vortex defect. Then, after some algebra, we obtain the result 

1 Yin + 1^2 

C{n2) = - E ^^n^[Mr) - i?n(r)], (11) 

^ n=0 
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1 — (x^ + + 2xr cos i 



a/2 



(10) 



where 



AJr) = i 



— 2r^ Inr, 

5 ' 



(4n+l) 2 _ 1 2n 
(n-l)(2n+3)' ' 



and 



— Inr, 



n = 1 



n = 



(12) 



(13) 



(2n+l)(4n+3) 2 _ (2n+l) 2n+2 „ n 
^ 2n(n+l)(2n+3)' 2n(n+l) ' ' 7^ U. 



This result is imphcit in an earher work of Bray and Humayun ||21| , who focused upon the 
singular part of this function. In the limit r — > 0, the singular terms in C(r) arise from 
Ai{r) and BqIj-), and can be computed as 

1 



Csing(ri2) = ^r^nr, 



(14) 



which gives rise to a power-law tail in the structure factor, S{k) ~ Ar^L'^ikL) ^, a result 



referred to as the "generalized Porod law" |22,15 . 



B. Case with /3 / 

We would like to undertake a similar asymptotic analysis in the general case with /5 7^ 0. 
As we are only interested in the limit r ^ 0, we again discard the step function on the RHS 
of Eq. d). In that case, we obtain 

^2^ 00 (^„T\n ,.l _ ^27r 

n\ Jo Jo 



C(ri2) 



IT 



n=0 



dxxe'"^^'' / d9{x + r cos 9) (x^ + + 2xr cos 6 



n-l 

2 



(15) 



We will separately consider the cases with n odd and n even, 
(i) n odd : 

We designate n = 2p + 1 and consider the angular integral on the RHS of Eq. (|l5l) : 
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I d9{x + r COS 9) {x'^ + r"^ + 2xr COS Oy 
Jo 



2xr^P I de(l + ^ + —coseY + 2r2p+i r dO cosO ( 1 + ^ + — cos o] " 



10 \ 
= 2xr^PIi + 2r2p+i/2. 



(16) 



The integrals Ji and I2 are obtained from Gradshteyn and Ryzhik [g^ , and the consohdated 
result is 



'2p+l 



[x, r) = 271 



^(P)2^2fc+1^2(p-fc)^^2 ^ ^P)^P+^)^2/c+1^2fc^^2^^2)p-2fc-l 



k=0 



k=0 



(17) 



where [y] refers to the integer part of y. The corresponding contribution to C{r 12) is 
Ciiru) = — - — (-1)^ ^, L dxxsm{qLx)In{x,r). 



n\ Jo 



(18) 



?i=l,3,5,.. 

The important feature here is that the above expression for Ci(ri2) only contains powers 
of r^. Therefore, the overall contribution to C(ri2) from this set of terms is analytic as 
r — i> 0. In the limiting case g = (/3 = 0), the above contribution is identically 0. 
(ii) n even : 

Next, let us consider the case with n even. We designate n = 2p, and the angular integral 
on the RHS of Eq. (|T3|) is 

f27r . _ _ \p-l/2 



hpix, ~ y (i6'(x + r cos 9) [x'^ + + 2xr cos 9 
Introduce p< = min(x, r) and p> = max(x, r) to obtain 



(19) 



hpix, r) = 2p?f-i / d9{x + r cos 0) 1 + ^ 



pI , 2p< 



p-1/2 



pI p> 



cos 9 



= 2p'P-\xh + rU). 



(20) 



The integrals and /4 can be computed in terms of hypergeometric functions as follows 
2|: 



/■? = -kF p, p: 1 

\^2 ^' 2 ^' p^ 



P< ) 

' -2 ' 



(21) 



and 
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TT 



1 



,2 

TT [p 



P 



P< 
P> 



2j p> I 2 



1 
2 
1 
2 



P> V2 



1 1 P< 



(22) 



We have simplified Eq. (p2D using the standard identity ||24 



(c — a — 1)-F(a, 6; c; + aF(a + 1, b; c; = (c — 1)-F(a, b;c — 1; z), 

with a = 1/2 -p, 6 = 1/2 -p, c = 2. 

Combining the expressions for I3 and /4, we obtain 



(23) 



hp{x,r) = 271 pI^-^ 



^F{^-p,^-p;l;^]+r^ {p-^)f{'^-P,^-P;2;^ 



,2 

pf/J 



(24) 



The corresponding terms in the correlation function are 

C2(ri2) = (1 - g') E (-1)"/'^^/ cixxcos(gLx) X 

n=0,2,4,.. ^' ^0 

„i /l — nl — n 



rpl~^P<{n - l)F I — , — ; 2; ^ 



V 2 



2 ' 'p|. 



;i-g2)(Ti + T2). 



(25) 



The singular contributions to C{r 12) as r — > arise entirely from C2(ri2), as Ci(ri2) is 
analytic in r. A considerable amount of algebra is involved in extracting the singular terms 
in Ti and T2. For the sake of brevity, we will only sketch the broad features of the calculation 
here. We have 



T^ = 2Y,{-IY 

p=0 



(2p)!r i-pj rn=0 



|2 



dxcos{qLx)x , (26) 





where we have used the standard expansion for the hypergeometric function The 
integral on the RHS of Eq. (p6D can be written as 
1 



h = rrr / dxcos(qLx)x 

j.2{m-p)+l Jq ' 



2m+2 _|_ ^2m 



dx cos{qLx)x 



-2(m-p)+l 



(27) 



The first term on the RHS of Eq. ( P7D is analytic as r ^ 0. The second term contributes 
singular terms only if m > p + 1 , yielding the result 
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frjLl2(m-p-i) 

Jg = (-i)^~P^^ In r + analytic terms. (28) 

[2(m — p — 

Replacing this in the expression for Ti, some algebra yields 

Ti = y y (-llP+^'+'-F^frT TT ^ r2(-+P+i) Inr + analytic terms. 

(29) 

A similar analysis for T2 yields 

T = V V r ^Y^m m'^'^-^ + (2^+1) ^^^^^^^^^ 

(2p)!(2m)!r(i_p)'(^ + p)!2(m + p+l) 

analytic terms. (30) 

We can combine the singular terms from Ti and T2 to obtain the singular part of C{r 12) 
as follows: 

aing(n2) = ^ E E (-^^" ^9 ^ 

2p=om=o (2p)!(2m)!r(i_p) (^ + ^ + 1)12 

(2m + 1) {2p + l)r2("^+P+i) In r. (31) 

We notice that the leading-order singularity is unchanged and continues to be Csing(ri2) — 
|r^lnr, as in the case with /3 = 0. However, there is now a sequence of sub-dominant 
singularities proportional to (gL)^r^lnr, (gL)^r^lnr, etc., and these become increasingly 
important as the length scale L increases. These sub-dominant terms in Csing(ri2) are 
reminiscent of the leading-order singularities in models with 0{n) symmetry, where n is 
Of course, in the context of 0{n) models, these singularities only arise for n < d 



even 



as there are no topological defects unless this condition is satisfied. In the present context, 
all these terms are already present for d = 2. The implication for the structure-factor tail 
is a sequence of power-law decays with S{k) ~ (qL)"^^"^^^^ L'^ / (kLY^"^"^, where m = 1,2, 
etc. Thus, though the true asymptotic behavior in c? = 2 is still the generalized Porod tail, 
S{k) ~ L'^{kL)~^, it may be difficult to disentangle this from other power-law decays. 
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The results presented in this section are of relevance in determining the small- distance 
behavior of the correlation function, or the large-wavevector behavior of the structure fac- 
tor. This is because small length-scales only probe individual defects. Nevertheless, as our 
numerical results in paper II will demonstrate, the single-spiral correlation function agrees 
with the correlation function for multi-spiral morphologies (obtained numerically) over a 
considerable range of distances. For even larger length-scales, we have to explicitly account 
for the modulation of the order parameter at defect-defect boundaries. We address this 
problem in the next section of this paper. 

III. UTILITY OF GAUSSIAN AUXILIARY FIELD ANSATZ FOR A 
MULTI-SPIRAL MORPHOLOGY 

The evolution in Figure 1 is characterized by a morphology with multiple spirals and anti- 
spirals. Initially, spirals and anti-spirals are attracted to each other and annihilate, thereby 
decreasing the defect density and increasing the inter-defect distance (or characteristic length 
scale). When the defect density is large, the spiral sizes are small and spirals are similar to 
vortices. Therefore, we expect an initial coarsening regime which is analogous to that for 
the XY model - both in terms of the domain growth law, L{t) ~ (t/lnt)^/^ and the 



morphology, as characterized by the correlation function ||Tl|J15| . This is in accordance with 
our numerical simulations, as we discuss in paper II. Distinctive effects of spirals are seen for 
length scales L > Lc, where Lc ~ - clearly, Lc ^ oo as g ^ (or j3 0). Furthermore, 
there is a repulsive spiral-antispiral potential beyond a certain distance, which prevents the 
annealing of all defects P,pO|. Thus, the evolving system "freezes" (in a statistical sense) 
into a multi-spiral morphology. This should be contrasted with the case of the dynamical 
XY model [a = j3 = in Eq. (|1|)), where we expect the zero-temperature system to continue 
coarsening as t ^ oo. 

A common theme in the characterization of dynamical evolution with a nonconserved 
order parameter is the introduction of a Gaussian auxiliary field (GAF) [^,^,[1^,^. Es- 
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sentially, the GAF ansatz takes the form V'(^)^) = F[m{f,t)], where the function F[m] is 
determined from the defect structure, and the complex field m (which measures the location 
relative to the defect core) is assumed to obey a Gaussian distribution. The zero-crossings 
of the field m correspond to the location of defect cores. The GAF ansatz enables a straight- 
forward computation of the correlation function for the field ip{f, t). However, the analytical 
justification for the GAF ansatz is meagre and its primary virtue appears to be that it works 
rather well in some situations [§]. 

Let us examine the utility of the GAF ansatz in the present context. The appropriate 
form of the ansatz for the CGL equation in the regime where the spiral structures are 
well-developed is (using Hagan's solution for the spiral defect) 

^(f,t) ~ ^ exi)[-i{ujt + q\m{f,t)\)], (32) 

'1 — -|- |m(r, t)p 



where u = (3{1 — g^); and we take {ipl ~ \m\ near the defect core {\m\ —>■ 0), in accordance 
with Hagan's solution. The field m(= mi + im2) is assumed to obey a Gaussian distribution 
with 

where cr^ = (mj(r, t)^); and the fields mi{f,t) and m2{r,t) are taken to be statistically 
independent of each other. 

Our numerical results show that the GAF ansatz in Eq. (^) is reasonable in the vicin- 
ity of defects. However, it is inappropriate for defect-defect boundaries, where the order- 
parameter amplitude \ip\ is often larger than y/1 — q'^. This is demonstrated in Figure 
4, which replots Figure 1 with defect locations marked by asterisks; and regions where 
\ip\ > \/l — q'^ marked in black. As discussed before, for early times (e.g., t = 25), the 
system evolution is governed by the interaction of vortices. Thus, the appropriate GAF 
ansatz should have maxd?/^!) = 1, as in the case of the XY model. For late times (e.g., 
t = 1000), the system has well-developed spirals. Nevertheless, the GAF ansatz for the 
order-parameter field is obviously inappropriate for large regions of space at these parame- 
ter values. For other values of (3, the same general arguments apply though there are changes 
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in the crossover time to spiral-mediated growth; and the fraction of spatial region where the 
GAF ansatz is unreasonable. 



Let us examine the validity of the GAF ansatz in regions where \tP\ < Vl ~ We can 
simplify the ansatz in Eq. ( |32D by defining the variable m' = me~^^, where (p = ut + q\m\. 
Then, we have the corresponding probability distribution for (say) m[ as 



/OO /"OO 
drrii / dm25{mi — m\ cos — sin <\))P (m\)P {mi) 
-OO J — OO 

]^ ^oo roo I 

— - / dm\ \ dm2S{m\ — rui cos (j) — irii sin (j)) exp - 

Tia^ J-oo J-oo \ 



ml + 



27ra2 J-oo ^ J-oo — " - — \ 2a2 

As usual, we transform {mi, mi) (|m|,^) to obtain 



(34) 



rfe5(m; - \m\ cos{9 - 0)). (35) 



Because of the periodicity of the function cos(6' — 0), the phase factor is inconsequential 
and 

and a similar distribution also applies for the variable m'g. 

Thus, we have the appropriate GAF ansatz (dropping primes) as follows: 

,^ . VI - q^m{f,t) 

Jl — q'' + |m(r, t)|"' 

where the variables mi{f,t) and m2{r,t) {m{f,t) = mi(r, t) +im2(r, t)) are Gaussian and 
independent of each other. The inverse relation between the variables ip and m is 

m(r,t) = , . (38) 

,Jl-q'-mr,t)\^ 

We want to examine the validity of the GAF ansatz numerically [0,0 in the context 
of the evolution depicted in Figure 1 (or Figure 4). The appropriate parameter values are 
/? = 1 and q ^ 0.306 [|1^] (see Figure 2). In Figure 5, we plot the single- variable distribution 
for the field mi{f,t), obtained directly from our simulation of the CGL equation using Eq. 
(^) in regions where \ip\ < Vl ^ Q^- The data in Figure 5 is obtained as an average over 5 
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independent runs for A^^-lattices, with = 512. (Details of our simulation will be provided 
in paper II of this exposition.) Figure 5(a) is a plot of P{mi) vs. nii from 4 different times 
- corresponding to the evolution pictures shown in Figure 1. In Figure 5(b), we have scaled 
variables and superposed the data for P(mi)cr vs. mi/ a, where a is obtained from the best- 
fit of the numerical data to the functional form in Eq. (^Bj). The data collapses onto a single 
master curve, which is reasonably approximated by the Gaussian form P{x) = 2 , 

denoted solid hue in Figure 5(b). 

Figure 5 has been obtained by focusing only on regions where < — g^, which is 
essentially equivalent to considering disjoint spirals, for which the correlation function has 
already been obtained in Section II. We have examined various ad-hoc methods of improving 



the GAF ansatz in Eq. (0). For example, one could set the saturation amplitude of the 
order parameter to its maximum value (iV'Isat — 1 for Figure 1), rather than {iplsat = ~ Q^- 
Figure 6 plots the resultant probability distributions, P{mi) vs. mi, with iV'Isat = 1- For 
early times (t = 25), the distribution has a Gaussian form, as expected from our analogy 
with the XY model. However, with the emergence of well-formed spirals, the distribution 
develops a double-peak and is clearly non-Gaussian. 

We have also studied some other possible ways of rectifying the GAF ansatz. We find 
that these ad-hoc approaches invariably result in non-Gaussian distributions for the auxiliary 
field. Perhaps a more honest approach should be based on the order-parameter field for 



spiral-spiral pairs |2^ as a function of two independent auxiliary fields - referring to distances 



from the centers of the two spirals. We are presently studying the utility of such an approach 
for characterizing the multi-spiral morphology. 

IV. SUMMARY AND DISCUSSION 

Let us conclude this paper with a brief summary and discussion of our results. We have 
undertaken a detailed analytical and numerical investigation of nonequilibrium dynamics in 
a special case of the complex Ginzburg-Landau (CGL) equation. Our results are described 
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in a two-stage exposition. This paper (referred to as I) constitutes the first stage of this 
exposition, and describes analytical results for the time-dependent correlation function. Our 
analytical arguments rely on the significance of spiral-defect structures in determining the 
morphology and evolution of the CGL equation from a random initial condition. 

In this paper, we describe results for the exact correlation function C(ri2) of a single spiral 
defect of size L, and undertake its asymptotic analysis in the limit ri2/L ^ but ri2/^ ^ 1, 
where ^ is the size of the defect core. We find that there is a sequence of singularities in 
this limit, which are reminiscent of singularities for defects with 0{n) symmetry, where n 
is even. However, the dominant singularity a.s r^/ L —>■ corresponds to the case of vortex 
defects, as expected. The implications for the large-wavevector tail of the structure factor 
are also discussed. 

We also investigate the validity of the Gaussian auxiliary field (GAF) ansatz in the 
context of multi-spiral morphologies. For early times (L < Lc ~ q~^), domain growth in 
the CGL equation is analogous to that for the XY model, whose domain growth law and 
correlation function are well understood [Q. For later times, we find that the simple GAF 
ansatz is not reasonable, as it is unable to account for order-parameter modulations in the 
defect-defect boundaries. We have attempted ad-hoc improvements of the GAF ansatz but 
these invariably result in non-Gaussian distributions for the corresponding auxiliary field. 
We are presently investigating the possibility of formulating a generalized GAF ansatz in 
terms of the order-parameter field for a spiral-spiral pair. 

More generally, the utility of the GAF ansatz arises from the summation over phases 
from many defects, which results in a near-Gaussian distribution for the auxiliary field. 
However, in the present context, the shocks between spirals effectively isolate one spiral 
region from the influence of other regions. As a matter of fact, the waves from other spirals 
decay exponentially through the shock and the phase of a point is always dominated by the 
nearest spiral. Therefore, we expect that the correlation function will be dominated by the 
single-spiral result - in accordance with our numerical results. 

The next paper in this exposition (referred to as II) will present detailed numerical results 
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for phase ordering dynamics in the CGL equation. In particular, we will focus upon the 
crossover from vortex-mediated dynamics (at early times) to spiral-mediated dynamics (at 
late times). Furthermore, we will compare our numerical results for the correlation function 
of the order-parameter field with the analytic form for a single-spiral defect presented in this 
paper. 

Before we conclude this paper, it is worth stressing that the results presented are easily 
adaptable to the general case of the CGL equation with a,/? 7^ 0. Again, the evolving 
morphology in a large region of parameter space is characterized by the presence and anni- 
hilation of spirals and anti-spirals 0. The results of the present paper apply directly in that 
case also, with minor modifications in the functional forms of the spiral solution in Section 
II. 
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FIGURE CAPTIONS 



Figure 1: Evolution of a small- amplitude random initial condition for the complex 
Ginzburg-Landau (CGL) equation with a = 0, P = 1. These evolution pictures were ob- 
tained from an isotropic Euler discretization of Eq. ([^), implemented on an A^^-lattice 
(A^ = 256) with periodic boundary conditions in both directions. The discretization 
mesh sizes were At = 0.01 and Ax = 1.0. The pictures show regions of constant phase 
9^ = tan~^(lmip/Reilj), measured in radians, with the following coding: 9^ G [1.85,2.15] 
(black); 9^, G [3.85,4.15] (red); 9^ G [5.85,6.15] (green). The snapshots are labeled by the 
appropriate evolution times. 

Figure 2: Plot of q{f3) vs. (3 for the 1-armed spiral solution of the CGL equation with 



a = 0. (cf. Figure 5 of Ref. 18 . 



Figure 3: Correlation function for the 1-armed spiral solution when /5 = 1 (g ~ 0.306). 
We plot C(ri2)/C(0) vs. r^/L for different spiral sizes, L = 10, 25, 50, 100 - denoted by the 
specified line-types. The results are obtained from a direct numerical integration of Eq. (^. 

Figure 4: Evolution shown in Figure 1, replotted to clarify the utility of the GAF ansatz 



in this context. The asterisks denote spiral centers; and regions where \^p\ > a/1 — are 
shaded black. 

Figure 5: (a) Plot of data for P[mi(r, t)] vs. mi(r, t) from 4 different times, t = 
25, 50, 100, 1000 - denoted by the symbols shown. The parameter values are identical to 



those in Figure 1. We use Eq. (38) to obtain data for mi(r, t) directly from the order- 
parameter field in our numerical solution of the CGL equation - considering only regions 
where {ipl < \/l — q'^. The data was obtained as an average over 5 independent runs for 
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Ar2-lattices {N = 512). 

(b) Scaled plot of data from Figure 5(a). We superpose data for P [mi{f,t)]a(t) vs. 
mi{f,t)/a{t), where a{t) is obtained from the best-fit of the numerical data to a Gaus- 
sian distribution. The solid line refers to the Gaussian function P(x) = -^e~^'^^'^. 

Figure 6: Plot of data for P[mi(r, t)] vs. mi(r, t) from times t = 25,50,100,1000 - 
denoted by the symbols shown. The parameter values and statistical details are identical 
to those for Figure 5(a). Data for mi{f,t) is obtained directly from the numerical data for 
7/'(r,t), using Eq. (^) with |?/'|sat = a/1 - q'^ replaced by |V'|sat = 1- 
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This figure "figl.gif" is available in "gif" format from: 
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